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ABSTRACT 

We have calculated the decay width of the process t — > gii including one loop QCD 
corrections. We found that the decay width of the such process could be comparable 
with that of the standard channel t — > bW'^, and, its QCD correction could enhance 
the widths over 30% in a very large mass range of the lighter stop. 



PACS number: 14.65.Ha, 14.80.Ly 



I. INTRODUCTION 



Top quark play an important role in the search for new physics due to its huge mass 
compared to other fermions in the standard model (SM). Besides its large Yukawa coupling 
to Higgs sector, top quark also anticipates strong interaction (QCD) with its supersymmetric 
partner through the top-stop-gluino vertex in a particle model of supersymmetry (SUSY). 
As a supersymmetric counterpart, a large mass splitting can emerge between the mass- 
eigenstates of the stops and could lead to a lighter physical stop, so that stop quark may be 
produced in our top quark factories. 

The experimental lower mass limit of stop is 86 GeV and the theoretical one could 
be more smaller 0. In most of conventional proposals stop will be produced in pair, then 
they need larger energy input in general. Comparing to stop pair production, one wonder 
whether stop could be produced associated with a light gluino in the decay of top quark. 
Such a channel can be accommodated in the single top events, and the top quark sector 
would become a nice place to study SUSY-QCD. 

The key point lies in that, whether the gaugino of SU{3)c is also sufficiently light. We 
would like to make an analyses on the situation of light gluino bellow. 

Historically, there had ever glimmered light gluinos in many scenarios with natural as- 
sumptions about the mechanics of SUSY breaking at higher energy scale. For example, the 
mass of gluino was resulted in 1 ~ 100 GeV typically if the SUSY breaking was transmit- 
ted to the observable sector by gravity (mSUDRA); and the masses of gluinos and squarks 
will arise smaller ( bellow GeV ) as a loop correction (next leading order effect) when that 
breaking is gauge-mediated ( GMSB ) ||^. Those light gluino scenario can explain sev- 
eral well-established phenomena, such as the anomalously slow running of SU{3)c coupling, 
anomalies in jet production as well as in the 7/(1410) 0. 

The difficulty is obviously that, present colliders do not favor light gluino. A lower mass 
bound 180 GeV has ever been published 0. Nevertheless, a controversial assumption that 
there is no gluino lighter than 5 GeV as well as the state-of-art usage of perturbative theory 
0, had made the constraint less convincing, especially leave a margin for a gluino light than 
5 GeV. The concomitant predictions of a light Higgs mass and a light chargino, have also 
been (at least marginally) ruled out at LEP II so has been the light gluino in mSUDRA 
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models. However, light gluino in GMSB models can survive in such experiment. Anyway 
all the direct searches for a light gluino have turned to a negative results since 1998. 

Then light gluino is cornered into bound states (for example ) to live if it is light 
enough. For a light gluino in a bound state, the well known KTeV experiment is an rather 
important one to give a constraint on its mass. KTeV result indeed has closed most of the 
place for a light gluino [|1^] and has given an impression that light gluino had been ruled 
out. Although, there are still two shreds that are not excluded by KTeV data in fact. The 
gluino could weight as ~ l.GGeV or as 5GeV [|^. Such light gluino can survive through 



another important E761 experiment too [11 



Although recent re-analyses |]T2| showed a consistency with SM and really eroded the 
attractiveness of light gluino in explaining slow running of SU{3)c coupling, such re-analyses 
relied too much on the measurement of low energy QCD, for example asiM^), which itself 
is difficult enough. 

Through above investigation, we have make it clear that, there is no counter- indication 
which can conclusively convince one to give up a light gluino ( few GeV or so). Recently, 
more new suggestion on light gluino within the framework of GMSB have been proposed 
r3| . Both scenario allow a gluino at few GeV and predicted the standard missing energy 



signature of SUSY may become not suitable. 

Then from the narrow window for light gluino, a noticeable partial width to stop and 
gluino might be detected in top quark decay when the statistics of the top events are im- 
proved. 

In the present paper, we will consider the decay process t — > gti including the QCD 



corrections, which possibility has been mentioned by Raby et. al JTSl- In experiments, the 
stop may be followed by the decays ti cXi or ti bxt [0], then its typical signal may 
be one hadronic jet plusing a large missing Pt- Instead of the detailed analysing the signal 
against possible background, we emphasize the possible anomalous branch ratio of top decay. 
In fact, in the SM BR{t bW~^) had been assumed as the dominant one for the observed 
top pair events at Tevatron ]l^. Thus, the non-SM branching ratio of the top quark which 
could be as large as the SM one is not excluded, and the searching for non-SM top decay 



are undergoing in both CDF and D0 |15]. 
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In the literatures, there have been several works [0 discussing the decays of stop or 
gluino into top quark, which have the same dynamics (interaction vertex) as the process 
discussed here, and these calculations offered a good references to check our calculations. 



II. CALCULATIONS 

The tree level diagram of the process t{p) ^(^1)^1(^2) is shown in Fig. 1 (a), and the 
decay width is 



3m' 



[m 



^ — m? + mf_^ — 2mgmi^ sin(26')), (1) 



where 



N = ^ {ml - {m-g + mi^y){m^ - {m-g - nitj'^), (2) 
and 6 is the stop mixing angle which is to be defined in next section. 

A. Virtual corrections 

The 0{as) virtual corrections arise from triangle as well as the gluino, quark and squark 
self energies diagrams as shown in Fig. 1 (6 ~ e) respectively. Through our calculations, we 
employed dimensional reduction regularization (DRR) to control all the ultraviolet diver- 
gences, which is widely adopted in the calculations of the radiative corrections in the minimal 
supersymmetric standard model (MSSM) since the conventional dimensional regularization 
violates supersymmetry with a mismatch between the numbers of degree of freedom of gauge 
bosons and gauginos [l^. At the same time, we adopted the on-mass-shell renormalization 
scheme [|l|,|20|. 

The one loop decay width can be expressed as 



-f2{{'mg-mt) -m^^){sme + cose)j , (3) 

where 
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/i = /r+/i', f2 = n+ft 



where /f are the contributions in the form of counterterms, 



f^ = V2{-hcos9-fnsm9), 

h = \zl + \zi + 5gjg, - tan + ^ 

fB^ = \z'^ + \zl + ^9s/gs + cot 656 + 



The renormahzation constants in Eq. ^ are, 

" 12^^? (sin0^(m| - m?)5o(0, m?, m|) + cos02(m| - m^)5o(0, m^, m|; 
— m^_Bo(0, A^, ml) + (^sin 6'^{m'^ — m^J — cos^^m^^ BQ{mf, rn^, mfj 

cos0^(m? — mf^) — sm6^m1^ BQ{ml,m^,ml^) 
+mf{—mg + Tn^^ — + 2mgmtsm[26))B{mf,mg,Tn^j 
+m1{—m?g + m| — — 2mgmt sm{26))B{ml, m?, m|) 

Zj, = Zi{6^n/2-6), 

/2| 2\D/2 2 2\ / 2| 2\D/2 2 2\ 

-(mg + m^JBoim-g, m^, m^J - {rrig + mgJBo{m~g, rrig, m^J 
■\-2m\{m\ - rr?g)B{m\, m\, m\^) + 2m\{m\^ - rr?g)B{m\, m\, m\^} 
+ (m| - mf)Bo{0,ml,ml) + (m| - mf)Bo{0,ml,ml) 

/2| 2 2\D/2 2 2\ / 2| 2 2\ d / 2 2 2 \ 

-{m-g + m^^ -m^)Bo{m~g,m^,m^J - [m-g + m^^ - mt)Bo{mg,mt,m^J 
+2m^K - -I - + 2-^-* sin(20))ij(m|, m?, m|) 
+2mg[m^^ — rrig — — 2mgmt sin(26'))i?(mg, , m^-^) 
-6m25o(0, 771% A^) + 24m|5(m|, , A^) 

% = ^ {-Bo{ml,ml, ml) + 5o(m|, A^ m| ) 

+ (m? - m| + mj^)5(m|, m?, ml) + 2m|5(m|, A^ m|) 
—2mgmt sm{26)B{ml_^ , m?, m^^ 
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We found the on mass shell definition for the counterterm of stops mixing, 



t2 11 



mtrrig cos(20) 



((So(m|, m?, m^2) - 5o(m|, m?, m^)) (7) 



is appropriate for our calculation, which has been adopted in |2T|]. One can exam that such 
a 66 is UV convergent, so that it is renormalization group ( scale ) independent at one loop 
level. As to the counterterm of the Yukawa coupling of SU{3)^, we prefer to the treatment 
in 
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^ y ^ {[A - log(4)]^ - i-Nc - — )}. (8) 

However there are fewer virtual particles decoupled in our scenario with light gluino and 
light squarks. 



= ^iM rl log + 1 log + 1 log + log ^ 

The ff in Eq. ^ are the contributions from the 3-points diagrams. 



fi = ""gg^^T^ - sin(2^))5o(m?,m|,m|) + sin(2^)5o(m?, m|, m|) 



+355o(m^, m^) + 18m^(3m^ + mt)^^^^ - 18{m'f - ml - ml)CP 
+ (-m^ - ml + m?)Ci^^ - 2m^t sm{26)C^^'^ + 2mt{-mg + mt + mt sm{26))Cf^ 
+ 18(-2A;i.p + 3m? + m^mt)^^ ^ + 18(-2A;i.p + 2m? + mymt)Cf^ - 2m^mtCf ^ 
—2mg{mg + {mg — mt) sin(26'))Ci'''' — 2mg{mt + {mt — m^) sin(2^^))Ci^'* 
+18(-4fci.p + 3m? + m'^t)cP + 18{-3ki.p + 2m? + m^)^^^) - 2{ki.p - m2)cf ^ 
— 2(m? — mgmt + {—ki.p + m? — mgmt + mf) sin(2^^))C2^'' 

— 2(— m^ + mgmt + (^i-P ~ ""^^ + fngf^t ~ ^1) sin (2^)) (72^"*^ , 

cos6' + sin6' / , . „ . 9 9 9 x , . 0/2 2 2 

]Bn(m. , m;;, m,- 



/2 = Qg^^2 (^(^ + sin (2^)) So (mi,m^,mj- J + sin(2^)z3ovm^ , m^, m^-^, 



-35^0(^0 A^ m^^) + 18mg(-3m^ + mt)Cj^^ + 18(m^ - m| - m?)Cj^^ 
-(-m^ - m| + m?)^^^^ - 2m(^ sin(2^)Ci^^ + 2mt{-mg ~mt + mt sin(2^))Cj^^ 
+18(2A;i.p - 3m? + m^mi)Cf ^ + 18(2A;i.p - 2m? + m^mt)Cf ^ - 2m3mtCf ^ 



+2mg{mg - (nig + mt) sm(26'))Cf ^ + 2mg{-mt + (nit + m-g) sin(26'))Cf ^ 
-18{-Aki.p + 3ml + m^t)Ci^^ - 18{-3ki.p + 2m? + m'^^)Ci^'^ + 2{ki.p - mf)C^^^ 
+2(m? + rrigmt + {ki.p — m? — rrigmt — mf) sin(2^))C2'^^ 

+2{—mf — rrigmt + {—ki-P + + mgmt + m^) sm{26))C2^^^ ■ (10) 

In the presentation of the formulas above, we have used notations, 

2 

- 7£; + /n(47r). 



ki-P 



mf — mf + 
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B{p^ml,ml) = dBo{p',ml,ml)/dp\ (11) 
pheavy -g decoupled virtual particles, which will not contribute to the scale evolution 



of as- The definitions of the scalar integrals Bs and Cs could be found in Ref. ||T9| , |22[| , and 
the indexes (1) — (5) of C functions have the variables 

/ 2 2 2 2 \2 2\/2 2 2 \2 2 2\/2 2 2 2 2 \2\ 

{m-,mt,m^^,mg,X , mj, (m^, , m^-^, A ,mg,m^^), {mg,mt,m^^,m^^,mt, X ), 
(m^, m^, m|, m^, m|, m?), (m?, m^, m|, m^, m|, m?), 

respectively. 

B. Real corrections 

The infrared divergence arise from the virtual massless gluon corrections are compensated 
by the real gluon bremsstrahlung corrections, i.e. the three-body decay 

t{p)^~g{k,Mk2)g{k). (12) 

The corresponding matrix element as given by the Feynman diagrams (Fig. 1 (f) ) is 

Mh = A^/2T:asU{kl){-ifabcTc^z^i{f + h+ mg){smePR - cos^P^) 
+TbTa-^{smePR - cos^Pl)(^i + ^2 + m*)/ 

+T,Tfe-i-(sin^P^ - cos^^Pi)(2A;2 + k).e}Uip), (13) 
2/C2.K 



where e denotes the polarization vector of gluon. Squaring the matrix element, performing 
the polarization and color sum over the square of the amplitude and integrating over the 
phase space yields the complete bremsstrahlung cross section as 



a: 



2A-nmt 



{48 \2ml{ml - m\ + - 2mgmt sin(2^))/ii + 2m?Ji + 



64 



-2mgmt sin(2^)/i] H 2m^ {—iTLg + m^^ — m^ + 2mgmt sin(2^))/( 

3 



00 



-2m^/o — Jq + 2mgmt sin(26') Jq 



64 



2m|(-m? + m| - 



+2mgmt sin(2^))/22 + 3m| - + 2mgmt sin(26'))/2 + / 



+ 3 



(2m^ — Arrigmi^ + 2m^-^ — 2m^ — Arrigmt sin(26')(m^ - m^-^ — m^))Ii^ 



02 



+ (m? + m| + m^)J2 + 2(m| - mfjl^ + 1 - 2m~gmt sin(2^)(/2 - Jq 
+24 -2(m? - m?^ + m^^)^ - 2mgmt sin(26')(m? - m| + m^))/oi 
+2(m^^ - m?)/o + 2(m?^ - + m-gmt sin(26'))(/i + Jq) - 21 



-24 



2(— m| + (m^^ — 2m'f^mf + ml + 2mgmt sin(26')(m? + m^^ — m^))/i2 



(14) 



Here the hard gluon has been included and the definition of function are 



jAB.. 

'-ab... 



■b^[p - fci - fc2 - k)- 



TT^ J 2k^2k^2k^ '^^ {±k.pa){±k.Pb)... 

where the minus corresponds top quark and plus corresponds gluino and stop|^. 



(15) 



III. NUMERICAL RESULTS 

In the MSSM the mass eigenstates qi and q2 of the squarks are related to the interaction 



eigenstates ql and g_R by |^| 



hA 


= W 




[<i2) 







with W 



cos 9q sin 9q 



\ 



V 



— sin 6q cos 6q j 



(16) 



* It should be noticed that (D.ll) and (D.12) in Ref. [^] must be corrected as following: the first 
term in the square bracket of (D.ll) must be replaced by the corresponding term in (D.12), and 

vice versa 
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Following the notation of , the mixing angle 6q and the masses m^^ ^ can be calculated 
by diagonizing the following mass matrices 



q 



=m}. + ml + ml cos 2/5 (/^^ - Cq sin^ Ou 



MIr 



rOfj ^ + + cos 2/5eg sin^ 9w 



(17) 



From Eqs. ^ and 17, and 9 can be derived as 



m 



_ 1 



tan 6' 



LL 



(18) 



For the numerical calculations, the renormalization scale is set at the mass of decay 
particle, fiR = rrit = 176.0 GeV as = 0.108. We also evaluated at fiR = mt/2as = .12 
and found the scale dependence very weak. The soft SUSY breaking mass terms are always 



m 



D 



m 



Q 



ms- For simplicity, we also define a variable 



naturally assumed as 
r = MlJMl^. 

In Fig. 1^, we show the branching ratios and the QCD relative corrections as function 
of m^^, assuming r = 0.9. The definition of branching ratio and relative correction are 

Br{t ^ gii) 



i+T(t^hw+)/v(t^gh) ^^"^ (5 = (Fi + r;,)/ro, respectively. /^From these curves, 
one can see that the branching ratios, which are not sensitive to the light gluino mass allowed 
by experiments, decrease with the increment of m^^ , which is the natural consequence of the 
shrink of the phase space. One also can see that, the decay width could be large enough to 
compare with the t hW^ channel in a large range of parameter space. On the other hand, 
the QCD corrections always enhanced the decay widths, for m^^ = 160 GeV, the corrections 
exceed 40%. 

In Fig. ^, we exam the result for its dependence upon the decoupled particles, by 
alternating the way for the mass of the lighter stop to vary. In these figures, r was scanned 
with fixed soft SUSY breaking mass term as 300 GeV. We can see these results are very 
close to Fig. 1^, which is due to the fact that the contributions from other squarks (relatively 
heavy) are less important. 
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In summary, we have calculated the decay width of the process t gti including the 
one loop QCD corrections. 

/^From our numerical examples, we can see that the QCD corrections could enhance the 
decay widths over 30% in a very large mass range of the lighter stop, and the decay widths 
of the process could be larger than that of the channel t —>■ bW~^. As a supplement of stop 
pair production at LEP 2000 or Tevatron, such a single stop ( top ) channel might make the 
top phenomenology more rich if it is allowed kinetic ally. 
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(a) 




FIG. 1. Feynman diagrams for the process t gti. 
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FIG. 2. Branching ratios and relative corrections as function of m^^, where r = 0.9. (a), the 
solid ( dashed) line represents the decay widths including QCD corrections for nig = 1.0 (5.0) GeV, 
and the dotted (dot-dashed) line denotes the tree level decay widths for nig = 1.0 (5.0) GeV; (b), 
the solid (dashed) line represent the relative corrections for nig = 1.0 (5.0) GeV 
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FIG. 3. see caption of Fig. but fixed the soft SUSY breaking mass term ms as 300 GeV. 
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